Ellipsoidal deformation of vertical quantum dots 
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Addition energy spectra at T of circular and ellipsoidally deformed few-electron vertical quan- 
tum dots are measured and compared to results of model calculations within spin-density functional 
theory. Because of the rotational symmetry of the lateral harmonic confining potential, circular dots 
show a pronounced shell structure. With the lifting of the single- particle level degeneracies, even 
a small deformation is found to radically alter the shell structure leading to significant modifica- 
tions in the addition energy spectra. Breaking the circular symmetry with deformation also induces 
changes in the total spin. This "piezo-magnetic" behavior of quantum dots is discussed, and the 
addition energies for a set of realistic deformation parameters are provided. For the case of the 
four-electron ground state at T, a spin-triplet to spin-singlet transition is predicted, i.e. Hund's 
first rule no longer applies. Application of a magnetic field parallel to the current confirms that this 
is the case, and also suggests that the anisotropy of an elliptical dot, in practice, may be higher 
than that suggested by the geometry of the device mesa in which the dot is located. 
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I. INTRODUCTION 

In atomic physics, the mean-field method describing 
the motion of electrons confined in the three-dimensional 
spherically symmetric Coulomb potential of the nucleus 
provides an impressively powerful tool to explain the 
chemical inertness and special stability of the noble gases. 
The well-known atomic shell structure is a consequence 
of the fact that the atomic levels Is, 2s, 2p, 3s, 3p, ... 
show a "bunchiness" in their distribution as a function 
of energy. Particular stability of the electronic system is 
reached when a bunch of such levels is fully occupied. If 
then one more electron is added, the electron configura- 
tion would involve a singularly occupied orbital from the 
next higher shell, and consequently, the system is then 
less stable. Shell filling is thus reflected by large maxima 
in the ionization energy for atomic numbers 2, 10, 18, 
corresponding to the nobel gas atoms He, Ne, Ar, ... . In 
the mid-shell regions, large level degeneracies occur as a 
consequence of the spherical symmetry of the confining 
potential of the atomic nucleus. The mid-shell levels are 
then filled according to Hund's rules, in particular max- 
imizing the total electron spin for half-filled orbitals |Q . 

A shell structure is not only unique to atoms, but ac- 
tually is a recurring property in finite fermion systems 
with high symmetry Q . It equally explains the occurence 
of "magic" proton and neutron numbers in the binding 
energies of nuclei, and more recently the discovery of 
"magic" atom numbers in metal clusters || - small ag- 



gregates of metal atoms in which delocalized valence elec- 
trons move in the positive charge background of the ions. 
Fundamentally, in contrast to atoms, however, both mid- 
shell nuclei and clusters deform their mean field rather 
than obey Hund's rules. 

The two-dimensional analogue to the atom with its 
static 1 /r radial Coulomb confinement due to the nucleus 
can be realized in small semiconductor devices: artifi- 
cial semiconductor atoms based on quantum dot technol- 
ogy. Clean, well-defined, and highly symmetric vertical 
quantum dots ( "islands" ) can now be made so small that 
the dot size is comparable to the Fermi- wavelength [|]-[| . 
Typical micrographs of micron-sized device mesas incor- 
porating these dots are shown in Fig. |l|. The lateral elec- 
trostatic confinement originates from side- wall depletion, 
and this (and the effective dot size) can be controlled or 
"squeezed" by the action of a Schottky gate wrapped 
around the mesa in the vicinity of the dot to the de- 
gree that the number of electrons trapped on the dot can 
be changed one-by-one. Also the few-electron regime is 
readily accessible, and then the resulting confinement is 
well approximated by a parabolic r 2 potential. Electron 
phenomena in related semiconductor quantum dot struc- 
tures continue to attract much attention |7|-^]. As they 
exhibit atomic-like properties, such as a shell structure 
and shell filling in accordance with Hund's first rule, the 
vertical quantum dots, whose heterostructure barriers are 
both abrupt and thin, can be regarded as artificial atoms 
whose ground and excited states can be probed electri- 
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cally by single electron tunneling spectroscopy in order 
to perform novel "atomic physics" experiments in the 
few-electron regime [p~0 1 1 . 

When an arbitrarily small bias, V, is applied across the 
dot between the metal contact on top of the device mesa 
and the substrate contact (these are often refered to as 
the source and drain contacts), the ground states of an 
N— electron quantum dot weakly coupled to the contacts 
can be investigated directly by monitoring the current 
flowing vertically through the dot at or below 0.3 K as 
the voltage on a single gate, V gi surrounding the dot is 
varied. When no current flows (Coulomb blockade), N is 
well defined. On the other hand, when current flows the 
number of electrons can oscillate between N and N + 1. 
With the gate, N can be increased one-by-one starting 
from zero by making V g more positive, so a series of sharp 
current peaks due to the charging of the dot (Coulomb 
oscillations) can be observed. For a large dot contain- 
ing many electrons, the Coulomb oscillations are usually 
periodic because the single electron charging energy is 
determined classically just by the total dot capacitance. 
For a dot containing just a few electrons both quantum 
effects reflecting the underlying symmetry of the con- 
fining potential, and the details of the electron-electron 
interactions become important as the dot size is reduced. 
This leads to modifications of the Coulomb oscillations, 
so they are no longer expected to be periodic pJlC|]. 

To date, we have mainly focused on the properties of 
dots in circular mesas which have diameters of typically 
0.4 to 0.7 microns. For a magnetic field parallel to the 
current, the measured ground states between T and 
about 4 T for N < 20 in these disk-shaped dot can be 
well accounted for by a single-particle picture based on 
the Darwin- Fock spectrum for a circular two-dimensional 
harmonic confining potential, a constant interaction, and 
corrections at T due to exchange, i.e. Hund's first rule 
H[To|. At higher fields beyond about 4 T, the evolution 
of ground states (and also the excited states) for N < 6 
can be understood in terms of many-body effects jll| . 

The main theme of this article concerns the effect of 
geometrically distorting a circular dot into an elliptical 
(anisotropic) dot. Previously, we have briefly reported 
some properties of elliptical dots ||[l2j- Here, we present 
a more detailed study of the addition energies and in- 
clude their magnetic-field dependencies. The experimen- 
tal data are compared to model calculations. We survey 
general trends, and examine basic assumptions about the 
nature of the deformed dots. 

A perfectly circular dot possesses full rotational sym- 
metry. This high symmetry leads to maximal level de- 
generacy of the single-particle two-dimensional states for 
parabolic confinement, and this emphasises atomic- like 
properties || . This level degeneracy at T for a circular 
dot is evident in the single-particle spectrum in Fig. |](a) , 
and consecutive filling of each set of degenerate states is 
directly responsible for the characteristic shell structure 
with "magic" numbers N =2, 6, 12, 20, ... . Further- 
more, Hund's first rule accounts for the parallel filling of 



electrons amongst half-filled degenerate states in a shell 
at numbers N —4, 9, 16, ... due to an exchange effect. 
Breaking the circular symmetry by deforming the lateral 
confining potential lifts the degeneracies of the single- 
particle levels present in a disk-shaped dot. This de- 
stroys the shell structure for a circle, and modifies other 
atomic-like properties [ fl3| ]. 

The sequence of spectra in Fig. || also introduces two 
key points in our subsequent arguments. Firstly, as the 
deformation is gradually increased, (a) to (d), degenera- 
cies of the single-particle states at T are generally re- 
moved. Nevertheless, accidental degeneracies can occur 
at certain "magic" deformations, e.g. (b) and (c), lead- 
ing to subshell closures, provided the confining potential 
is still perfectly parabolic. The resulting patterns, how- 
ever, are very different from that for the circular case, 
(a), and in practice may be hard to observe. Secondly, 
a weak magnetic field parallel to the current can also 
induce level degeneracies in both circular and elliptical 
dots when single-particle levels cross at finite field, but 
here too, any shell structure at a particular field is of a 
lower order and less apparent than that for the circle at 

While illustrative, ultimately any modelling of the be- 
havior of real dots must go beyond a system of N non- 
interacting electrons confined by a two-dimensional har- 
monic oscillator, i.e. a single-particle picture, as em- 
ployed to generate the spectra in Fig. || Q, and include 
Coulomb interactions which can lift certain degeneracies 
at T. Numerical diagonalization of the full Hamiltonian 
matrix has recently been successfully employed to calcu- 
late basic electronic prope rties of dots with anisotropic 
confining potentials |14|,[l5| . Such "exact" numerical cal- 
culations, however, are limited to only a few confined 
particles. In order to study dots confining a larger num- 
ber of electrons we apply spin-density functional theory 
at T. This powerful technique, which explicitly incorpo- 
rates the electron-spin interactions, has lead to a number 
of interesting predictions for the ground state structure of 
quantum dots, although there is a continuing discussion 
as to the interpretation of so-called spin-density waves 
(SDW) fl6j-|2l|l . Both "exact" numerical calculations and 
spin-density functional theory predict subtle changes in 
the addition energy spectra, and transitions in the spin- 
states as deformation is varied - even for a weak defor- 
mation. An example of the latter is the breakdown of 
the conditions for which Hund's first rule applies for four 
electrons, and this marks a transition from a spin-triplet 
to a spin-singlet configuration, i.e. states are consecu- 
tively filled by spin-up and spin-down electrons. 



II. EXPERIMENTAL SETUP 

The vertical quantum dots under focus in the fol- 
lowing are fabricated by electron-beam lithography, 
and a two step etching technique to make circu- 
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lar or rectangular sub-micron mesas from one special 
GaAs/Al .22Gao.78As/In .o5Gao.95As/Alo.22 Ga .7sAs/- 
GaAs double barrier heterostructure (DBH) . Full details 
of the device fabrication, and the material parameters 
are given elsewhere PH^,|lO|l- A single Schottky gate is 
placed around the side of the mesa close to the DBH. We 
discuss one circular mesa with a nominal top contact di- 
ameter, D, of 0.5 /Ltm (W), and three rectangular mesas 
with a top contact area (L x S) 0.55 x0.4 /xm 2 (X), 0.65 
xO.45 /im 2 (Y), and 0.6 x0.4 /mi 2 (Z). L(S) is the nom- 
inal dimension of longest (shortest) side of the top con- 
tact. Fig. ^] shows typical scanning electron micrographs 
of a circular mesa, and a rectangular mesa taken imme- 
diately after the depositon of the Schottky gate metal 
surrounding the mesa. For the rectangular mesas, an 
intuitively simple way to classify them is to define a geo- 
metric parameter, f3, to be the ratio L/S. For X, Y, and 
Z respectively (3 is nominally 1.375, 1.44 and 1.5. Due 
to a slight isotropic undercut resulting from the light wet 
etch during the formation of the mesa [|||, the area of 
the mesas, as revealed by the micrographs, is a little less 
than that of the top contact, so realistic values for j3 are 
estimated to be about 5% larger than the values quoted. 

Fig. [I] also schematically shows the slabs of semicon- 
ductor between the two Alo.22Gao.7sAs tunneling barri- 
ers, and the resulting dots bounded by the shaded deple- 
tion region for the circular and rectangular mesas. The 
thickness of the Ino.05Gao.95As slab is determined by the 
separation between the well defined heterostructure tun- 
neling barriers (approximtely 100A). The slab is suffi- 
ciently thin that all electrons are in the lowest state in 
the vertical direction parallel to the current. The lateral 
confining potential due to the side wall depletion further 
restricts electrons to the center of the slab, thus defining 
the dot region. We note that in our devices, the extent 
of the lateral depletion region in the vicinity of the dot is 
largely determined by the electron density in the n-doped 
GaAs regions above and below. 

The lateral harmonic confining potential of the dot in 
the circular mesa has circular symmetry of a sufficiently 
high degree that degenerate sets of states can systemat- 
ically form in the disk-shaped dot l|. These states can 
be labelled by the quantum numbers (n,l), where n is 
the radial quantum number (=0, 1, 2, ...), and I is the 
angular momentum quantum number (= 0, ±1, ±2, ...). 
Each state can hold a spin-up electron and a spin-down 
electron. At T the 2n + |2| + 1-th shell is made up of 
2n + \l\ + 1 degenerate single-particle states. Each de- 
generate set of states can be regarded as a shell of an 
artificial atom, and this is the origin of the 2, 6, 12, 20, 
... "magic" numbers. The first shell consists of the (0,0) 
level, the second shell of the (0,1) and (0,-1) levels, the 
third shell of the (0,2), (1,0) and (0,-2) levels, and so 
on. For the circular dots we typically study, the lateral 
electrostatic confinement energy seperating these degen- 
erate sets of single-particle states, Eq, can be as large as 
5 meV in the few-electron limit pT| . Neglecting an ar- 
bitary constant, the energy of single-particle state (n,l) 



is (2n + \l\ + 1)Eq. The effective lateral diameter can 
be "squeezed" from a few thousand Angstroms for N of 
approximately 100 down to A for N = by making 
the gate voltage more negative ||,|l(|[ll| . We stress that 
crucially the "squeezing" action of the gate, and indeed 
application of a magnetic field parallel to the current, 
preserves the circular symmetry of a disk-shaped dot. 
Consequently, atomic-like properties should be particu- 
larly robust and evident in circular dots. 

For a rectangular mesa, the lateral confining poten- 
tial of the dot is expected to be elliptical-like due to 
rounding at the corners provided the number of elec- 
trons in the dot is not too large (in which case it may 
be more rectangular- like with rounded corners), or too 
small. Right at "pinch-off", (N — > 0), it may even 
become more circular-like, i.e. the elliptical-shape of 
the confining potential may be changing in a complex 
way [|]|lj]- Assuming the confining potential is per- 
fectly parabolic, we can choose to characterize the "ellip- 
ticity" by a deformation parameter, 5 — E$/El. Here, 
Es(El) is the confinement energy at T along the mi- 
nor (major) axis (E$ > El). The states in the el- 
liptical dot are now labelled by the quantum numbers 
(n-L,ns), where Ul (ns) is a quantum number (=0, 1, 
2, ...) associated with the energy parabola along the 
major (minor) axis Again neglecting an arbitary 

constant, the energy of single-particle state (riL,ns) is 
(n L +l/2)E L +(n s +l/2)E s . 

For a perfectly circular mesa, we can trivially general- 
ize our definition of the deformation parameter so that 
6 = j3 = 1. On the other hand, for the rectangular mesas, 
there is no simple correspondance between /3, a ratio of 
lengths characteristic of the top metal contact which is 
independent of gate voltage (or N), and S, a ratio of en- 
ergies characteristic of the dot in the mesa which is in 
fact dependent on the gate voltage (or equivalently N), 
i.e. "accidental" degeneracies at "magic" deformations 
will be hard to see over an extended range of N, and 
in any case may be lifted if the confinement potential is 
not completely parabolic. Nevertheless, at this stage, we 
start by assuming that (3 is a measure of <5, and thus one 
might expect Sz > Sy > Sx > Sw- We are not saying 
that 6 = p for the ellipses, and indeed even for the sim- 
plest possible model of uniform depletion spreading due 
to the action of the gate, we would expect (3 to underes- 
timate 8. We furthermore assume in the following model 
calculations, for simplicity, that the "squeezing" action 
of the gate does not alter 6. We will examine these as- 
sumptions in light of the experimental and theoretical 
data presented. Note that the application of a magnetic 
field parallel to the current effectively reduces 5 as seen 
by the confined electrons in the limit of a very high field, 
where it approaches unity. 
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III. ADDITION ENERGY SPECTRA FOR 
CIRCULAR AND DEFORMED DOTS 



In Fig. ||, the change (formally the second difference) in 
the electro-chemical potential, fi(N+l)—fj,(N) — A2(AQ, 
which can also be regarded as a capacitive energy [fL6| , 
is plotted as a function of electron number, TV, up to 
N = 17 for (a) W, (b) X, (c) Y, and (d) Z at T. The 
traces are offset vertically by 3 meV for clarity. The 
Nth current peak position in gate voltage, V s , at a very 
small bias (<C 1 mV), i.e. measured in the linear con- 
ductance regime, at or below 0.3 K reflects n(N), the 
electro-chemical potential of the ground state for N elec- 
trons, or equivalently the "addition energy" to place an 
extra electron on a dot with N — 1 electrons. A2(iV) then 
mirrors directly the spacing in gate voltage between the 
N + 1th and the iVth current peaks 0. A 2 (N) is ac- 
tually the half-width of the Nth Coulomb diamond, the 
diamond-shaped region in the V — V g plane in which cur- 
rent is blocked between the Nth and the N + 1th current 
peaks. A2 contains contributions from the single-electron 
charging energy and changes in the single-particle energy, 

e q i& 

At T, for the circle W, i±2(N) is strongly depen- 
dent on N , and a very clear characteristic shell struc- 
ture is evident in Fig. |](a) |10|. Particularly large peaks 
(N = 2, 6, 12), and relatively large peaks (N = 4, 9, 16) 
are indicated. The result from a local spin density ap- 
proximation (LSDA) calculation discussed below is also 
included for comparison ^0|. N =2,6, and 12 are the first 
three "magic" numbers for a circular two-dimensional 
harmonic potential which mark completion of the first 
three shells (containing respectively 1, 2 and 3 degener- 
ate zero-dimensional single- particle states or equivalently 
2, 4 and 6 electrons). The peaks at N = 4, 9, 16 arise as 
a consequence of exchange effects which are enhanced at 
half-full shell filling with same-spin electrons for the 2nd, 
3rd, and 4th shells respectively (To). This shell structure 
should be clear (and this is generally the case in practice 
for N < 20) as long as: i) the two-dimensional lateral po- 
tential remains radially parabolic, and rotationally sym- 
metric to a fairly high degree, ii) Eq is comparable to, 
or larger than, the Coulomb interaction energy, and Hi) 
the effect of screening is not significant. 

For the circular mesa W, it is also evident that as 
N is decreased, A 2 (iV) generally becomes larger due to 
the increase of the Coulomb interaction when the dot 
is "squeezed". This observation also holds for the rect- 
angular mesas, but there are no prominant maxima at 
A2(2,6,12). The shell structure for the disk-shaped dot 
has now become disrupted or "smeared out" , and this can 
be attributed directly to the lifting of the degeneracies of 
the single-particle states that are present in a circular 
dot ]q jl0| , ^2| |. In other words, deformation kills the shell 
structure for a circle, and even quite a small deformation 
can make a big difference. This is evident from the three 
traces, (b) to (d) in Fig. 0, but there are major difficulties 



in discussing specific details. As noted earlier, in prac- 
tice, right at "pinch-off", 5 may actually tend towards 
unity §,|l^|, but more generally 8 « f3 may be unreliable. 
Also, even for two circular dots which have a clear shell 
structure in the few-electron limit, the absolute values of 
A 2 (N) can vary from dot to dot, i.e. the precise details 
are device dependent, and beyond the third shell only 
a few devices show the expected behavior clearly [[10| . 
Lastly, even if S could be determined accurately, A 2 (N) 
strictly speaking can only be fairly compared if the "ar- 
eas" of the dots are comparable, as in the classical limit 
A 2 (N) is determined by the overall dot capacitance [||. 
Based on the nominal sizes of the mesas, and in line with 
the trends of the "pinch-off" gate voltage as identified 
by the position of the first current peak, elliptical dots 
X, Y, and Z respectively may have "areas" 1.1, 1.5 and 
1.2 larger than that of the circular dot W. Thus to sen- 
sibly discuss details, like the spin-states, even generally, 
we first calculate A2(N) at T for a range of 6 values 
in line with those suggested by the (3 values of the mesas 
X, Y, and Z. We then compare with, and look for pat- 
terns in, the experimental data at T, before looking at 
the magnetic field dependence for confirmation of trends, 
and whether 6 ~ (3 is reasonable. 



IV. MEAN-FIELD MODEL FOR CIRCULAR AND 
ELLIPTICAL QUANTUM DOTS 

We next aim to model the changes due to the defor- 
mation of the lateral confinement to the shell structure 
of the quantum dots at T by applying the methods of 
spin-density functional theory (SDFT). We will briefly 
address different aspects of the spin structure relevant to 
the deformed quantum dots. 



A. The method 

To obtain the ground-state energies and densities for 
N electrons confined in an externally imposed potential, 
we solve the spin-dependent single-particle Kohn-Sham 
(KS) equations [§§ 



2m* 



(1) 



in a plane- wave basis to avoid any symmetry restrictions. 
In Eq.(|j]), the index a accounts for the spin (| or J,), and 
r = [x, y). The effective mean-field potential, V° s (r), 
contains contributions from the external harmonic con- 
fining potential, the Hartree potential of the electrons, 
and the functional derivative of the local exchange- 
correlation energy, for which we use the approximation 
of Tanatar and Ceperley |23| (see also 17 2C|] for details). 
The electrostatic confinement due to the lateral depletion 
region imposed by the side wall and the Schottky gate is 
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approximated by a two-dimensional anisotropic harmonic 
oscillator with frequencies u> x — and uj v = uj/y/S, 

V ext (x,y) = ^m*Lj 2 (sx 2 + ±y 2 ^j . (2) 

The ratio of the oscillator frequencies, S = lo x /ui v , thus 
defines the ratio of semiaxes of the ellipsoidal equipo- 
tentials. We impose the constraint, which is 

equivalent to conserving the area of the quantum dot 
with deformation ppj . The x and y-axes are indicated 
in the schematic diagram for the elliptical dot in Fig. [I]. 
With this convention, the above defined Eg and El re- 
spectively correspond to Huj x and htu y . In the model we 
present the dot is assumed to be well isolated from its sur- 
roundings, so any effects due to the presence of the gate 
and the neighboring conducting regions are neglected. 
Likewise, screening and non- parabolicity effects inside 
the dot, which become more important for large N, are 
not considered. 

For 5 = 1, a circular shape for the quantum dot is 
obtained, whereas 5 > 1 corresponds to an ellipsoidally 
deformed quantum dot. The strength, to, of the exter- 
nal parabolic confinement leading to an average particle 
density, no = l/(iTr 2 ), in a circular dot is approximated 
by uj 2 — e 2 1 (4neoem*r3yf]\[) [^7). Minimizing the energy 
density functional by self-consistcntly solving the above 
KS equations, Eq. (|l|), ground state energies, E(N, 6), are 
obtained for different electron numbers and deformation 
parameters. Full technical details are given elsewhere 
0,^O|, and here we report only the results. We em- 
phasise that from recent measurements, it is clear that 
as N increases the confinement weakens in such a way 
that the particle density tends to a constant . This is 
implicit in our model, as for any given value of r s , the os- 
cillator frequency to, and the related frequencies u x and 
w v , decrease with increasing N . S is also kept constant 
for simplicity, although S is expected to vary with N in 
practice. 

Although strictly speaking the dot is located in 
Ino.05Gao.95As, we take for values of the effective mass, 
m* , and dielectric constant, e, those for GaAs - namely 
0.067 and 13.1 respectively. There are no fitting param- 
eters in the equations, and only a suitable choice for r s 
is required to generate the addition energy spectra. The 
value of r s = l.5a B used in the model calculations is real- 
istic as the value estimated experimentally for a circular 
quantum dot is 1.3 to 1.4a* B Q. a* B = h 2 (47re e) / m* e 2 
is an effective atomic unit, which for GaAs is about 103A. 
r s = l.ba* B in the model presented here corresponds to 
an effective confinement energy, Eq, for N = 1 of about 
5.7 meV. This value is consistent with the upper limit of 
Eq observed in practice (about 5 meV), and justifies the 
Eq=3 meV value as a reasonable average for calculating 
the simple single-particle spectra shown in Fig. ^ for the 
first ten levels. 

We point out that the SDFT calculations described 
here, as well as those performed by Hirose and 



Wingreen ]2l| , are strictly two-dimensional, so the 
strength of the Coulomb interactions may be overesti- 
mated, i.e. the possibility of charge spreading out in both 
the x-y-planc, and along the vertical direction parallel to 
the current to minimize the Coulomb energy is neglected. 
Equivalently, anisotropic extension of the electron wave- 
functions along the major axis is ignored Q. In prac- 
tice, screening by the metal contacts surrounding a dot 
is also believed to reduce the influence of Coulomb inter- 
actions. The three-dimensional model of Lee et al. |l(| 
does incorporate self-consistent solution of the Poisson 
equation into a SDFT calculation, but because they use 
different expressions for the exchange-correlation energy, 
and considerably higher values for Eq, El, and E$ (up 
to 20 meV), a direct comparision with their results is 
not easy. Nevertheless, they find that electrons strongly 
confined in the vertical direction have a very strong two- 
dimensional character, and both approaches lead to the 
same qualitative conclusions. Namely, the distinct shell- 
structure for a circle, as well as the spin-states, at T 
are strongly modified with deformation. 



B. A 2 (iV) for elliptical dots: Results from LSDA 
calculations 

We now make a simple comparison between the exper- 
imentally measured traces for the change in the electro- 
chemical potential, A2(A), with those modelled theo- 
retically. Fig. [| shows A2(iV), derived from the self- 
consistent ground-state energies, E(N,S). The ener- 
gies are obtained by self-consistently solving the KS- 
equations starting from different initial guesses for the 
effective KS-potential for the spin-up and spin-down par- 
ticles. The initial potentials are chosen completely arbi- 
trarily by just putting small random numbers on to the 
lattice points. The calculations are started from four such 
guesses. For two of them, the spin-up and spin-down ini- 
tial guesses are shifted in value in order to to search for 
states with non-zero total spin for even-A. This is impor- 
tant in order to find the ground state amongst all possible 
spin configurations with a high degree of certainty p"7| , p0"[ . 

The lowest trace in Fig. ^ gives A2(N) for the circu- 
lar dot (6=1, i.e. zero deformation). As expected, the 
circular- shaped confinement produces a spectrum with 
the familiar shell structure for a two-dimensional har- 
monic oscillator, with shell closures at the "magic" num- 
bers 2, 6, 12, and 20. At the average particle density cor- 
responding to r s — l.5a B , these "magic" numbers arise 
from large gaps at the Fermi surface and paired spins 
in each non-degenerate level, so the total spin is zero 
(5 = 0). We note that within this mean-field model, spin- 
density wave (SDW) states are not expected for these 
particular spin-zero states p7|^0| . 

In Fig. ||(a), for 5 = 1, the experimental and theoret- 
ical traces can be directly compared. The agreement is 
strikingly good, given that no parameters are fitted to 
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reproduce the experimental data. Not only are the prin- 
cipal peaks 2, 6 and 12 well reproduced, but the rela- 
tively large peaks at 4, 9, and 16 for the high-spin states 
at half-shell filling are also clear jl|,|§. For N = 4, 
Hund's first rule correctly predicts the calculated 5=1 
spin-triplet state in which spins are aligned in the two 
highest partially occupied degenerate single-particle lev- 
els (n,l)=(0,l) and (0,-1), rather than the 5 = spin- 
singlet state in which the paired-spin electrons reside in 
either the (0,1) or (0,-1) levels. 

Deforming the confinement slightly by changing the de- 
formation parameter to 5 = 1.1 (see trace (b) in Fig. |^), 
the calculation still predicts fairly clear shell closures at 
N = 2,6 and 12. These numbers can still be consid- 
ered as "magic", but the actual values of A 2 (2, 6, 12) are 
noticeably suppressed, because degeneracies have been 
lifted J(| . The N — 20 peak has become very weak. Also 
values of A 2 (7V) neighboring N =2, 6, and 12 start to be- 
come comparable to the values for N = 2, 6, and 12, i.e. 
there is less contrast. Overall, the shell structure is much 
less pronounced compared to that for the circle. Already 
it is clear that even a very small deviation from perfect 
circular symmetry can have a very noticeable effect even 
when single-particle level degeneracies are lifted by just 
a small amount. 

As the deformation increases further, the pronounced 
peaks for N =2, 6, 12, and 20 evident for the disk-shaped 
dot are further suppressed. This is a simple consequence 
of the removal of the level "bunching" with deformation. 
Even for the cases where "accidental" subshell closures 
occur at certain "magic" deformations (e.g. S = 1.5 and 2 
as seen in Fig. ||) , the reduced separation between degen- 
erate single-particle energy levels (El) would make any 
shell structure less clear to observe, and the sequence 
of "magic" numbers would be very different (e.g. for 
6 = 2 it would be 2, 4, 8, 12, 18, ...) compared to those 
for 6 = 1. From Fig. [| we can see essentially that for 
5 > 1.2, the circular shell structure has been completely 
eliminated. Traces (a) to (f) thus illustrate the dramatic 
destruction of the familiar shell structure for a circular 
dot with deformation. 

Also apparent is that a systematic one-to-one corre- 
spondance of A 2 (iV) between traces (b) to (d) in Fig. || 
and traces (b) to (f) in Fig. ^ is impossible to make. Al- 
though the experimental data for mesa X partly resem- 
bles the theoretical data for 6= 1.1 to 1.3, the data for 
mesas Y and Z do not seem to resemble that for S > 1.3, 
except perhaps for a weak tendency to oscillate between 
even-iV and odd-./V. We have already stated many rea- 
sons why, in comparision to a circular dot, a good corre- 
spondance between experiment and theory for the ellip- 
tical dots is less likely. We stress that ultimately, except 
for circle W, 5 is not known, and equating 5 with (3 may 
not be reliable. To progress we must look for other clues. 

Theoretically, Fig. [| shows that there are transitions 
in the ground state spin-configurations with deforma- 
tion p0[ . The total spin, 5, is identified by different 
symbols in the figure. These transitions are particularly 



numerous for, but are not restricted to, the even-iV sys- 
tems, and are clearly very sensitive to the actual value of 
the deformation. For example, in the case of N = 6 elec- 
trons, the total spin is predicted to change from 5 = 
(i.e. a paramagnetic state) at 5 = 1, through an 5 = 
SDW state, to S = 1 at 5 = 1.5 — an indication of 
"piezo-magnetic" behavior |2^j2^], i.e. changes of the dot 
magnetization with deformation. Although experimen- 
tally we are not in a position to differentiate between an 
5 = "normal" state and an 5 = SDW state showing 
a spatial variation in the polarization as a consequence of 
broken spin symmetry in the internal coordinates p6|]- 
indeed the interpretation of a SDW is still debated in 
the literature |21]- the SDFT calculations described here 
predict that the latter becomes more prevalent for even- 
N systems as S increases, particularly for small average 
particle densities [ ^7|j20| . 

Another interesting, and in practice the simplest inci- 
dence we can focus on, is what happens to the N = 4 
ground state. The inset in Fig. |I] shows A 2 (Af = 4) ver- 
sus deformation up to 5 = 1.5. Starting with the circular 
dot, Hund's first rule gives a total spin of 5 = 1 for the 
triplet state favoring spin alignment of the two electrons 
in the second shell rather than a total spin of 5 = for 
the singlet state in which the spins are paired. As the 
deformation is initially increased, the energy separation 
between the two levels (nL,ns)= (1,0) and (0,1) - the 
two originally degenerate levels (n,l)= (0,1) and (0,-1) in 
the second shell of the circular dot - increases (see (a) 
and (b) in Fig. ||), and so the spin-triplet state becomes 
progressively less favorable. A 2 (4) continously decreases 
with 5, and at a value between 1.2 and 1.3, a spin-zero 
state (actually predicted by the SDFT described here to 
be a SDW) appears, i.e. a spin triplet-singlet transition 
is expected. For higher values of 5 beyond this transition, 
A 2 (7V = 4) starts to increase. 

Other recent calculations employing numerical diago- 
nalization for elliptical dots moderately deformed up to 
5 — 2 have also predicted that A 2 (TV) is sensitive to defor- 
mation, and that the spin-states can be modified 0,E5[. 
Those calculations, for N up to 10, and performed at 
T with Eq=3 meV, also reveal a spin triplet-singlet 
transition at 6 « 1.2 for N = 4, and, more generally, 
a consecutive filling of states by spin-up and spin-down 
electrons at higher deformation is favored. 

Inspection of Fig. ^ gives values of A 2 (7V = 4) for 
mesas W, X, Y, and Z respectively of 3.1, 2.7, 3.1, and 
2.5 meV. Whilst it is reassuring that these energies lie 
in the range predicted by SDFT, it is tempting to at- 
tribute, for a 5 value equated to the (3 value, the appar- 
ently anomalously low value for mesa Z to sample spe- 
cific fluctuations, and say that the trend for mesas W, 
X, and Y is consistent with that predicted in Fig. [| (in- 
set), i.e. TV = 4 is a spin- triplet for W, and a spin-singlet 
for X, Y, and Z. However, as we do not really know 5 
for the elliptical dots, we can not even be confident that 
the actual 5 values lie in the 5—1.0 to 1.5 range, i.e. the 
values might be higher, or even that the order X, Y, and 
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Z for increasing deformation as suggested by the j3 values 
is correct. Fortunately, we can apply a i?-field, and as we 
will shortly show this goes a long way to resolving these 
difficult issues. 

In case the actual 8 values for the elliptical dots ex- 
ceed 1.5, traces (g) and (h) in Fig. [| respectively show 
A2(iV) for the higher deformation parameters 6 = 2 and 
<5 = 3.2. We have no reason to believe that (5 experimen- 
tally will be exactly 2 or exactly 3.2, but the numbers are 
representative of the two situations where, respectively, 
many or no single-particle levels are degenerate at T for 
non-interacting electrons, as illustrated by the spectra in 
Fig. ||. As expected, traces (g) and (h) show no circular- 
like shell structure, and no particularly large values of 
A2(-/V). Indeed, apart from the "classical" background 
trend, i.e. A2(A) increasing as A decreases, there is lit- 
tle one can say about the traces except for A > 5 there is 
a tendency for a weak even-odd oscillation in A2 (A) , and 
this oscillation is perhaps clearer for larger S. The model 
here actually predicts small peaks for odd-A, and small 
valleys for even-A. For odd-A the spin-state is nearly 
always S=l/2, and for even-A the spin-state is usually 
5 = (SDW). At least for 6 = 2, where in the single- 
particle picture there can be accidental degeneracies at 
T (see Fig. || trace (c)), one might naively expect some 
non-zero even-A spin-states, but it is possible that in the 
model calculations, for the parameters given, the interac- 
tions modify the spectrum so dramatically that expected 
degeneracies are lifted reducing the visibility of any po- 
tential shell structure, e.g. A = 6 and 10 are predicted 
here to be S = (SDW) rather than S = 1 as might 
be expected from Hund's first rule. On the other hand, 
for 6 = 3.2, where in the single-particle picture there are 
no accidental degeneracies at T (see Fig. || trace (d)), 
perhaps surprisingly some non-zero even-A spin-states, 
for example for A=12 and 16, are predicted - this too 
may be due to interactions. 

For Y and Z, the A2(A) traces in Fig. |^ seem to show 
a weak tendency to oscillate between a slightly larger 
even-A value, and a slightly smaller odd-A value, and 
this oscillation seems clearer for Y than for Z. For the 
moment we do not try to account for the clarity of this 
oscillation in dots Y and Z, but try to explain the origin 
of the oscillation, although we are now being forced to en- 
tertain the idea that S for Y and Z may be much larger 
than 1.5. Starting from the over simple single-particle 
picture with a fixed confinement energy, and then includ- 
ing a constant interaction which is the same for even-A 
and odd-A, a larger even-A value is expected because 
only A2 (even-A) can contain a finite contribution due 
to the single-particle energy level spacing. A slightly 
more advanced model, which is more realistic in prin- 
ciple, would be to have a constant interaction for odd-A 
(next electron added to an S = 1/2 state already con- 
taining one electron) that is stronger than the constant 
interaction for even-A (next electron added to an empty 
state). If the former is larger than the latter plus the 
single-particle spacing (more likely in practice as A in- 



creases), a weak tendency to oscillate between smaller 
even-A and larger odd-A could occur. This pattern is 
what the SDFT calculations predict in Fig. 4 for 5 = 2 
and S = 3.2. The fact that A 2 (A) for Y and Z is often a 
little larger for even-A than odd-A should not be taken 
to mean that the constant interaction model is more ac- 
curate. Rather the Coulomb interactions may not be so 
strong in practice, due to screening by the leads for ex- 
ample, as those in our model- a model that also does not 
include the self-consistent calculation of the electrostatic 
confining potential. Indeed, in the SDFT calculations 
of Lee et al. [[l6| , the electrostatic confining potential is 
much stronger (e.g. Es=20 meV, El=10 meV), and they 
find that A2(A) is generally a little larger for even-A 
than for odd-A. Finally, we note that eventually, for 
a much stronger deformation (e.g. 6 exceeding 10), the 
addition energy spectrum would become smoother as it 
tends towards that for a quasi-one-dimcnsional quantum 
wire @|^. 



V. MAGNETIC FIELD DEPENDENCE 

Application of a magnetic field is a powerful tool with 
which to identify the quantum numbers of states in our 
vertical quantum dots ||[l0]jll]]. Fig. || instructively 
shows the expected evolution of the first ten single- 
particle energy levels with £?-field up to 6 T for a cir- 
cular dot (<5=1), and for elliptical dots with 6 = 1.5,2, 
and 3.2. The energy level spectra are calculated accord- 
ing to the simple single-particle model employed by Mad- 
hav and Chakraborty jl3| in which Coulomb interactions 
are neglected, and the confining potential is assumed to 
be perfectly parabolic. The spectrum for the circular 
dot is the familiar Darwin-Fock spectrum for a circu- 
lar two-dimensional harmonic confining potential. The 
confinement energy for the circular dot, Eq, is taken 
to be 3 meV, in practice a reasonable average value in 
the few-electron limit, and is assumed to be indepen- 
dent of A. The confinement energies for the elliptical 
dots are simply derived from the relation E^Es=EqEq. 
For the case of the circle and the 5 = 3.2 ellipse, quan- 
tum numbers (n,l) and (n_L,ns) respectively for some of 
the states we discuss are indicated. Each single-particle 
energy level can accomodate a spin-up and spin-down 
electron, so current peaks should normally come in pairs 
in a constant interaction model neglecting exchange [fl3[ . 
"Wiggles" in the position of pairs of current peaks are 
expected because the _B-ficld induces crossings between 
single-particle states |||ll],[l3|] . The first lowest energy 
"wiggle" originates from the crossing marked by a black 
triangle in each of the four spectra. A2 (even-A) is ex- 
pected to be strongly dependent on B- field as it can con- 
tain contributions from single-particle energy level spac- 
ings, whereas A2 (odd-A) is essentially independent of 
_B-field at weak-field, and is determined only by the ef- 
fect of Coulomb repulsion. Any detailed discussion on the 
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actual B-field dependence of the current peaks requires 
the inclusion of Coulomb interactions jllj . The four cal- 
culated spectra nevertheless clearly serve to demonstate 
three simple points: i) the S-ficld lifts all degeneracies 
present at T at the "magic" deformations, e.g. 5 =1, 
1.5, 2, ... (<5=3.2 is not a "magic" deformation); ii) a B- 
field can always induce degeneracies at finite field when 
single-particle levels cross, provided the confinement po- 
tential is perfectly parabolic; and in) as S increases, the 
single-particle energy level spacing generally decreases 
(<El)_ 

Fig. H shows the £?-field dependence, for a weak field 
applied parallel to the current, of the Coulomb oscilla- 
tion peak positions for the circular mesa W, (a), and the 
rectangular mesas X, Y and Z, (b) to (d). The data con- 
sists of current vs. V g traces taken at a very small bias 
(<C 1 mV) at different B-fields at or below 0.3 K. 

For circle W, only the third, fourth, fifth and sixth 
current peaks (belonging to the second shell at T) are 
shown. The pairing of the third peak with the fifth peak, 
and the fourth peak with the sixth peak from T to 0.4 
T, as opposed to the more usual pairing of the third peak 
with the fourth, and the fifth peak with the sixth (due 
to consecutive filling of electrons into spin-degenerate 
single-particle states) for B>0.4 T, is a consequence of 
Hund's first rule: the N = 4 state is a spin-triplet so two 
parallel-spin electrons fill the two different but originally 
degenerate states (n,l) = (0,l) and (0,-1) in the half-filled 
second shell §,[l0). For B>0.4 T, the fifth and sixth 
peaks, as a pair, first move up, as indicated by the thick 
arrow, and then start to move down at about 1.4 T due to 
the crossing of the single-particle states (n,l)= (0,-1) and 
(0,2). This lowest single-particle level crossing, which is 
also clear in Fig. ||(a), is marked by a black triangle. The 
spins of the added electrons are also shown pictorially at 
T and 2 T. 

To explain why Hund's first rule is obeyed in a sim- 
ple way, we can introduce an energy, Eex, to repre- 
sent the reduction in energy due to exchange between 
electrons in the half-filled second shell, and this is es- 
timated to be about 0.7 meV for circle W §,0. The 
N = 4 triplet-state is thus lower in energy than the 
N = 4 singlet-singlet state by Eex, and as a conse- 
quence A 2 (3),A 2 (5) < A 2 (4) by about 2E E x- This 
exchange-related effect persists in a weak B-field as long 
as the splitting between states (0,1) and (0,-1) is less 
than Eex- At 1.4 T this splitting exceeds Eex, and the 
ground state becomes a spin-singlet, i.e. there is a B-ficld 
induced triplet-singlet transition. 

For rectangles X, Y, and Z, the first ten current peaks 
are shown in Fig. |],(b) to (d). Peaks are paired, and there 
are no obvious deviations close to T for N — 4 which 
can be attributed to exchange effects, i.e. Hund's first 
rule. Quantum numbers (nL,ns) of the single-particle 
states are assigned, and the first up-moving pair of peaks 
is marked by a thick arrow. With increasing deformation, 
the first up-moving pair of peaks, and the lowest energy 
single-particle level crossing- identified by a black trian- 



gle in each of the Fig. || spectra- are simply expected 
to move systematically to higher N (or equivalently to 
higher energy) 

For the elliptical dots, normal peak pairing, even from 
T, occurs so Hund's first rule is not obeyed. This 
suggests that the spin-state for N = 4 is a singlet. 
The exchange effect is maximal for a circular dot at 
N — 4 because the (n,Z)=(0,l) and (0,-1) states are de- 
generate, but with deformation these states become the 
(riL,ns)=(l,0) and (0,1) states in an elliptical dot which 
are split at T. This energy splitting, 7, increases with 
S. If 7 < Eex at T, exchange can still operate to lower 
the energy, and thus the N = 4 ground state remains 
a spin-triplet. On the other hand, if 7 > Eex at T, 
the energy gain due to exchange is not sufficiently large 
to compensate for the spliting, so normal pairing occurs. 
Thus, as 5 increases, we can expect a triplet-singlet tran- 
sition at some critical deformation ]l2] j. Note that Eex 
itself decreases with increasing deformation, as it has its 
maximum value only when the orbitals involved have the 
same symmetry. This transition is clear in the inset of 
Fig. ^ according to SDFT, and has also been predicted 
by exact numerical diagonalization [ fL^JTB] . The tell-tail 
pattern in the trend of A 2 (A^ = 4) at T with deforma- 
tion should be an initial decrease while the state remains 
a spin-triplet, a turning point at the transition, and a rise 
thereafter when the state is a spin-singlet. As noted be- 
fore, it is hard to judge from the absolute values of A 2 (4) 
at T alone shown in Fig. || whether the N — 4 state is a 
triplet or singlet. A 2 (4) can be relatively large either side 
of the turning point if either Eex or 7 is large, i.e. a large 
A 2 (4) can mean Hund's first rule is operating for nearly 
degenerate states, or there is a large separation between 
non-degenerate states. This potential ambiguity is ap- 
parent when we see that A 2 (4) for circle W and ellipse Y 
are essentially equal, so it is vitally important to examine 
carefully the B- field dependence. The absense of devi- 
ations to the normal peak pairing at N = 4 in Fig. [jj, 
traces (b) to (d), nevertheless does apparently confirm 
that S is indeed greater than 1.2-1.3 which is in line with 
the j3 values for mesas X, Y, and Z. For completeness, 
we note that normally we probe the spin-states in our 
high symmetry dot structures via the orbital effect with 
the B-field parallel to the current. The spin-states in 
the elliptical dot X have also been confirmed directly by 
measuring the Zeeman effect alone by applying a _B-field 
perpendicular to the current, and the results are again 
consistent with a spin-singlet interpretation for N = 4 

& 

The next most striking feature about traces (b) to (d) 
in Fig. |^ is the position of the first up-moving pair of 
peaks. For mesas X, Y, and Z respectively, it is the 3rd, 
5th, and 4th pair of peaks. As revealed by the sequence 
of spectra in Fig. |^, in a simple single-particle picture, the 
first up-moving state is (n£,ns)=(0,l), which is actu ally 
the lowest energy state of the second Landau-level [ ]13| . 
Inspection of these calculated spectra shows that this 
state is, in the weak- field limit, from the bottom, the 3rd, 
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4th, and 5th state respectively for 1 < 6 < 2, 2 < S < 3, 
3 < S < 4. Thus, starting from no deformation, the first 
up-moving pair of peaks should go from the 3rd to 4th, 
4th to the 5th, ... at certain "magic" deformations as 
6 is increased. Remembering that Coulomb effects are 
neglected in this simple picture, and that in practice 5 is 
expected to vary with N, nonetheless, with these simple 
arguements it looks as if 1 < S < 2 for X, 3 < 6 < 4 
for Y, 2 < 5 < 3 for Z. If we believe this, then even 
though ellipses X, Y, and Z are all deformed beyond the 
triplet-singlet transition, we are forced to conclude the 
following: i) S can be much higher than that suggested 
by the (3 values (especially for Y and Z); and ii) the or- 
dering given by increasing (3 values may not reflect the 
true ordering in S, i.e. the deformation in Y seems to 
be stronger than in Z, so the true sequence may be W- 
X-Z-Y for the four mesas considered. Given our earlier 
comments, the former is not so unexpected since we have 
no independent way of measuring S, but the latter is per- 
haps more surprising. 

If the true ordering of the ellipses is X, Z, and Y, even 
though the reason for the deformation in Y being stronger 
than in Z is unclear to us, at least other attributes of Y 
and Z, and trends in Fig. g and Fig. || are consistent with 
this interpretation. For instance, A2(Y = 4) for a W-X- 
Z-Y ordering respectively of 3.1, 2.7, 2.5, and 3.1 meV is 
more in line with the predicted trend shown in the inset 
in Fig. ||, although the value for Z still seems a little low. 
This reordering does not contradict our earlier conclusion 
that, with the absense of deviations to normal pairing, 
N = 4 is a spin-singlet state for all three ellipses. Thus 
A2(4) increases after the triplet-singlet transition, be- 
cause with deformation the degeneracy of single-particle 
states is strongly lifted. The reversed ordering of ellipses 
Y and Z, as well as higher 8 values than suggested by 
the j3 values, also fits with the observations made earlier 
about traces (c) and (d) in Fig. ||, and traces (g) and 
(h) in Fig. ^. Both mesas show a tendency for A2(N) to 
oscillate between slightly higher and slightly lower values 
respectively for even-iV and odd-iV and this seems more 
pronounced for Y than Z. 



VI. CONCLUSIONS 

We have experimentally and theoretically investigated 
the effect of ellipsoidal deformation on the shell structure, 
addition energies, and spin-states in vertical quantum 
dot atoms on going from circular- to rectangular-shaped 
mesas. The familiar and distinctive shell structure as de- 
termined from the addition energy spectra at T for the 
circular dot is absent in the elliptical dots, and even small 
deviations breaking circular symmetry have a dramatic 
effect. Measurements with a magnetic field applied paral- 
lel to current confirm that the N — 4 spin state at T has 
undergone a transition due to the moderate deformation: 
for the circular dot it is a spin-triplet in accordance with 



Hund's first rule when the second shell is half-filled, and 
for the elliptical dots it is a spin-singlet. These observa- 
tions are in agreement with recent theory, as well demon- 
strated here by the application of spin-density functional 
theory at T with a wide range of deformation parame- 
ters. The S-field dependence strongly suggests that the 
anisotropy of an elliptical dot in practice can be signifi- 
cantly higher than that given by simply considering the 
geometry of the mesa in which the dot is situated. In 
the future it will be interesting to experiment with even 
more strongly- and extremely-deformed dots (to clarify 
for instance the existance of S = SDW states), quasi- 
one-dimensional wire-like dots |2(],^7j , and possibly other 
exotically-shaped dots, for example ring-shaped dots p7[ , 
and triangular-shaped dots 0,15). Finally, with these 
goals in mind, ultimately better control and in-situ ma- 
nipulation of the lateral potential geometry of a quantum 
dot is highly desirable, and this may be achieved by fully 
exploiting a multiple-gated vertical single electron tran- 
sistor we have recently developed P§,E9 . 
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FIG. 1. Typical scanning electron micrographs show a cir- 
cular mesa with top contact diameter D, and a rectangular 
mesa with top contact area L x S (L > S) taken imme- 
diately after the deposition of the Schottky gate metal sur- 
rounding the mesa. The slight undercut due to the light wet 
etch is clearly visible. Schematic diagrams depict the slabs 
of semiconductor between the two tunneling barriers, and the 
resulting circular and elliptical shaped dots bounded by the 
shaded depletion region for the circular and rectangular mesas 
respectively. The current flows vertically through the dots in 
the direction indicated. For the rectangular mesa, the energy 
parabola along the major and minor axes are included and 
the respective confinement energies, El and Es(El < Es), 
are marked. 

FIG. 2. Magnetic (B-) field dependence up to 6 T of the 
first ten single- particle energy levels (each level can hold a 
spin-up and spin-down electron) for a circular dot with 5 = 1, 
(a), and for elliptical dots with 5 =1.5, 2, and 3.2, (b) to (d). 
The energy level spectra are calculated with the simple for- 
malism employed by Madhav and Chakraborty, as explained 
in the text. The spectrum for the circular dot is in fact the 
familiar Darwin-Fock spectrum for a circular two-dimensional 
perfectly harmonic confining potential. The confinement en- 
ergy for the circular dot, Eq, is taken to be 3 meV, a reason- 
able average value in the few-electron limit, and is assumed to 
be independent of N. The confinement energies for the ellipti- 
cal dots are simply derived from the relation ElEs = EqEq. 
For the circular dot and the 5 = 3.2 ellipse, quantum numbers 
(n, I) and (riL,ns) respectively, for some states are given. At 
T accidental degeneracies are also evident for the 5 = 1.5 
and 5 = 2 ellipses, but not for the 5 = 3.2 ellipse. The black 
triangles mark the position of the first (lowest energy) sin- 
gle-particle crossings 



FIG. 3. Change in the electro-chemical potential, A2(N), 
as a function of electron number, N, up to 17 at T for cir- 
cular mesa W of diameter 0.5 fi, (a), and three rectangular 
mesas X, Y, and Z of area 0.55 x 0.4 /im 2 , 0.65 x 0.45 /jm 2 , 
and 0.6 x 0.4 fj,m 2 respectively, (b) to (d). For W, 5 = 1, and 
for X, Y, and Z respectively, /3 is nominally 1.375, 1.44, and 
1.5. The traces are offset vertically by 3 meV for clarity. For 
circular dot W, a clear shell structure is observed. Peaks due 
to full (N = 2,6,12) and half-full (N = 4,9,16) shell filling are 
numbered, and the fit given by local spin- density approxi- 
mation (LSDA) is included, as discussed in the text and also 
shown in Fig. 
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FIG. 4. Model calculations for the change in elec- 
tro-chemical potential, A2(N) within spin-density functional 
theory. The different traces correspond to zero, weak and 
moderate deformation parameters 8 — 1.0 to 8 = 1.5, (a) to 

(f) , and higher deformation parameters 8 = 2 and 8 = 3.2, (g) 
and (h). The traces are offset vertically by 1 meV for clarity, 
and there is an additinal 1 meV offset between traces (f) and 

(g) . Traces (a) to (f) illustrate well the dramatic destruction 
with deformation of the familiar shell structure for a circular 
dot. The total spin, S, for different deformations and elec- 
tron numbers are identified by different symbols as defined 
in the figure. Note that with increasing deformation, 5 = 
spin-density wave (SDW) states are predicted to become more 
prevalent for even-JV. The inset shows A2(JV = 4) versus 8. 
For 8 < 1.2 and 8 > 1.3 respectively the N = 4 ground state is 
expected to be a spin-triplet (S=l) and spin-singlet (S = 0). 



FIG. 5. Magnetic (B-) field dependence of the Coulomb os- 
cillation peak positions for mesas W, X, Y and Z, (a) to (d). 
The B-field is parallel to the current. The data consists of 
current versus V g traces for different B-fields which are offset, 
and rotated by 90 degrees. A less negative gate voltage cor- 
responds to higher energy. For circular dot W, only the 3rd, 
4th, 5th, and 6th current peaks belonging to the second shell 
are shown. The unusual pairing of the 3rd peak with the 5th, 
and the 4th peak with the 6th from to 0.4 T as opposed to 
the more usual pairing for B>0.4 T is evident, and is related 
to Hund's first rule. For B>0.4 T, the 5th and 6th peaks, as 
a pair, first move up, as indicated by the thick arrow, and 
then start to move down at about 1.4 T due to the crossing 
between the single-particle states (n,Z)=(0,-l) and (0,2). The 
spins of the added electrons are shown pictorially at and 2 
T. For elliptical dots X, Y, Z, the first ten current peaks are 
shown. Peaks are paired, and there are no obvious deviations 
close to T for N — 4 which can be attributed to exchange ef- 
fects. Quantum numbers (riL,ns) of the single-particle states 
are included, and the first up-moving pair of peaks is marked 
by a thick arrow. As deformation increases, the single-particle 
spectra in Fig. 2 show that this up-moving pair should move 
to higher N. 
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